Cosmological and Astrophysical Parameter Measurements with 21-cm Anisotropics 

During the Era of Reionization 
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We study the prospects for extracting cosmological and astrophysical parameters from the low 
radio frequency 21-cm background due to the spin-flip transition of neutral Hydrogen during and 
prior to the reionization of the Universe. We make use of the angular power spectrum of 21-cm 
anisotropics, which exists due to inhomogeneities in the neutral Hydrogen density field, the gas 
temperature field, the gas velocity field, and the spatial distribution of the Lyman-a intensity field 
associated with first luminous sources that emit UV photons. We extract parameters that describe 
both the underlying mass power spectrum and the global cosmology, as well as a set of simplified 
astrophysical parameters that connect fluctuations in the dark matter to those that govern 21-cm 
fluctuations. We also marginalize over a model for the foregrounds at low radio frequencies. In 
this general description, we find large degeneracies between cosmological parameters and the astro- 
physical parameters, though such degeneracies are reduced when strong assumptions are made with 
respect to the spin temperature relative to the CMB temperature or when complicated sources of 
anisotropy in the brightness temperature are ignored. Some of the degeneracies between cosmolog- 
ical and astrophysical parameters are broken when 21-cm anisotropy measurements are combined 
with information from the CMB, such as the temperature and the polarization measurements with 
Planck. While the overall improvement on the cosmological parameter estimates is not significant 
when measurements from first-generation interferometers are combined with Planck, such a combi- 
nation can measure astrophysical parameters such as the ionization fraction in several redshift bins 
with reasonable accuracy. 

PACS numbers: 98.80.Es,95.85.Nv,98.35.Ce,98.70.Vc 



I. INTRODUCTION 

The 21-cm spin-flip transition of neutral Hydrogen, ei- 
ther in the form of an absorption or an emission rela- 
tive to Cosmic Microwave Background (CMB) blackbody 
spectrum, provides one of the best ways to study the in- 
tergalactic medium during and prior to reionization 
With frequency selection for observations, the 21-cm line, 
in principle, provides three-dimensional tomography of 
the reionization era as well as a probe to the dark ages 
where no luminous sources are present after recombina- 
tion 0, S H i, !, • The exact physics associated with 
the reionization process is still largely unknown, though 
it is strongly believed that UV photons from first lumi- 
nous sources are responsible for it Q. 

The 21-cm background could reveal some details of the 
reionization process including when Lyman-a photons 
first began to appear, through the Wouthuysen-Field ef- 
fect M , when the gas was heated by an X-ray radiation 
field [10l | . and the formation of first virialized halos or 
"mini halos" [ll|. The 21-cm background also captures 
the underlying mass fluctuations, since inhomogeneities 
in the gas density field are expected to trace those of the 
cold dark matter (CDM), and, thus, fundamental cos- 
mological parameters that define the linear density field 
power spectrum [(| [l2l [l3j . If imaged with adequate 
resolution, in principle, 21-cm anisotropies can be used 
to "de-lens" CMB B-mode polarization maps and to im- 
prove the energy scale of inflation well below the limits 



allowed with lensing analysis in CMB data alone [1J|. A 
major concern for all these studies is the extent to which 
dominant foreground signals can be separated with the 
21-cm signal from reionization extracted out; techniques 
have been proposed and developed based on the smooth- 
ness in the frequency space of foreground signals while 
21-cm itself varies rapidly [la. Ilq. Il7l Il8l|. 

Two recent studies considered the possibility to mea- 
sure cosmological parameters by making certain simpli- 
fying assumptions related to the sources of fluctuations 
in the 21-cm brightness temperature such that the 21-cm 
aniso trop ies are only related to those of the gas density 
field jl2l. Il3j. In a more general scenario, 21-cm bright- 
ness temperature fluctuations are sourced by a variety of 
inhomogeneities including the gas density, the gas tem- 
perature, the gas velocity, and the Lyman-a intensity 
field [!, El H3l The assumption of a single source of in- 
homogeneity involving the gas density requires the sce- 
nario that the spin temperature of neutral Hydrogen is 
significantly higher than CMB temperature, such as due 
to heating by an X-ray radiation field [10J . While such an 
assumption has been generally used in the description of 
the 21-cm background fluctuations P. fl6{. it is not clear 
exactly when the gas heated, and the time scale related to 
the heating process, or even if the gas heating precedes or 
follows the complete reionization. Over a wider range of 
redshift, it is likely that the intergalactic medium (IGM) 
is in a state where spin temperature is not vastly different 
from that of the CMB and that Lyman-a photons from 
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first sources contribute to 21-cm anisotropics as well [20|. 

Here, we study anisotropics in the 21-cm background 
during an intermediate regime where Lyman-a photons 
exist and perturbations are introduced to the 21-cm 
brightness temperature. We focus on the measurement 
of cosmological parameters in such a scenario, by tak- 
ing into account additional sources of anisotropy. Since 
processes that govern 21-cm fluctuations are now domi- 
nated not just by the inhomogeneities in the gas density 
field, but rather by a large number of sources, we also pa- 
rameterize a set of "astrophysical quantities" that relate 
the underlying physics to fluctuations in the 21-cm sig- 
nal. These astrophysical parameters capture (in various 
combinations) the mean gas temperature, the ionization 
fraction of the IGM, and the scale-dependent bias fac- 
tors that relate fluctuations in the density, temperature, 
and the Lyman-a radiation intensity field to that of the 
underlying CDM power spectrum. 

Here, we discuss the difficulty to measure cosmolog- 
ical parameters in the presence of uncertainties associ- 
ated with astrophysical parameters. We also consider the 
combination of 21-cm information with those from the 
CMB, especially using Planck temperature and polariza- 
tion observations. The combination of 21-cm anisotropics 
and CMB data breaks certain parameter degeneracies, 
but this is not adequate to significantly improve cos- 
mological parameter measurements over Planck alone. 
The combination does improve, however, parameter es- 
timates from the astrophysical side that are degenerate 
with cosmological quantities. We suggest that, in the 
limit where 21-cm anisotropics need to be studied in a 
more generalized context than suggested in the litera- 
ture under a narrow scenario, upcoming first-generation 
low radio frequency interferometers, such as the Mileura 
Wide-field Array (MWA 1 ) and the Low Frequency Array 
(LOFAR 2 ), are, at best, suitable to study astrophysics 
during reionization; these interferometers may allow the 
determination of a single parameter involving the prod- 
uct of the neutral fraction and the neutral gas bias factor 
at the level of 10% in several redshift bins prior to com- 
plete reionization. 

With significant improvements on the instrumentation, 
such as with the Square Kilometer Array (SKA 3 ), it could 
be that one can make reasonable improvements in cos- 
mological parameter measurements, especially parame- 
ters such as the tilt of the primordial power spectrum or 
the cosmological constant, beyond the level one can reach 
with CMB measurements from Planck. In general, how- 
ever, given the parameters involved from the astrophysics 
side, and their degeneracy with cosmological parameters, 
it is unlikely that 21-cm observations alone will be com- 
petitive with the precision one can reach with CMB. The 
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mean signal of the 21-cm background also captures as- 
trophysics related to the reionization process, and, if in- 
dependently measured, could further help with breaking 
of degeneracies in the astrophysical parameters; we do 
not include such information here as planned interferom- 
eters, by design, are not sensitive to the mean background 

dUll. 

In addition to discussing a more general scenario for 
21-cm background anisotropies, our analysis also consid- 
ers complications associated with the analysis and inter- 
pretation of the 21-cm signal. Following Ref. [r|, we 
include a model for the residual foregrounds after fore- 
ground removal based on the multifrequency technique 
given the smoothness of the contaminants in frequency 
space @. We also include cross-correlation of 21-cm sig- 
nals between frequency bins, which cannot be ignored due 
to fluctuations generated by the gas velocity field that 
has a larger correlation length than density fluctuations. 
This cross-correlation of the signal between channels re- 
duces the efficiency to which foregrounds can be cleaned 
and the cosmological information captured by 21-cm fluc- 
tuations relative to the case where bins are considered to 
be independent. 

The paper is organized as follows: in Qll Al we present 
a general outline of the 21-cm signal from various redshift 
ranges since the recombination. In ijllBl we discuss a gen- 
eral formulation of the 21-cm signal anisotropy and fre- 
quency correlations, with the focus on the angular power 
spectrum of 21-cm fluctuations measured in terms of the 
brightness temperature relative to the CMB tempera- 
ture. In i]IIII we discuss two model descriptions of the 21- 
cm fluctuations as relevant for this calculation, focusing 
on the appearance of Lyman-a photons that couple gas 
temperature to that of the spin temperature of neutral 
Hydrogen, and an era when the Universe has undergone 
some reionization with gas temperature heated signifi- 
cantly above that of the CMB; the latter is the scenario 
generally studied in the literature since it both provides a 
simple description of the 21-cm anisotropies and is at the 
low redshift ranges targeted by first generation interfero- 
metric observations. In QIV\ we discuss the experimental 
setup assumed in our analysis, noise contributions, as 
well as the foreground model at low radio frequencies. 
As part of our model fitting procedure, we marginalize 
over all parameters related to the foregrounds at the same 
time as we extract cosmological and astrophysical infor- 
mation. In SjV]we present results related to cosmological 
and astrophysical parameter measurements using 21-cm 
anisotropy measurements and conclude with a summary 
of our results in § VI. Throughout the paper, we make use 
of the WMAP-favored ACDM cosmological model [H[. 



II. 21-CM BACKGROUND: MEAN SIGNAL 
AND THE ANISOTROPIES 

In this section, we will first discuss the mean brightness 
temperature of the 21-cm signal relative to that of the 
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blackbody CMB and then the fluctuations in the bright- 
ness temperature due to inhomogcneities of the density, 
temperature, and velocity fields, among others. We di- 
vide our discussion of the 21-cm signal into several red- 
shift ranges between us and the last scattering, depending 
on the dominant physics of interest and the reionization 
history of the Universe. 



A. Brightness Temperature 

When traveling through a patch of neutral hydrogen, 
the intensity of the CMB radiation will change due to ab- 
sorption and emission. The corresponding change in the 
brightness temperature, T 2 i, as compared to the CMB at 
an observed frequency v in the direction n is then 



T 2 i(n, v) 



l + z 



(1) 



where T$ is the temperature of the source (the spin tem- 
perature of the IGM), z is the redshift corresponding 
to the frequency of observation (1 + z = 1/21/ V, with 
V21 = 1420 MHz) and T 7 = 2.73(1 + z)K is the CMB 
temperature at redshift z. The optical depth, r, of this 
patch in the hyperfine transition [2J| is given in the limit 
of k B T s » hv 2 i by 



3c 3 HAio 7lHI 



lGkv^Ts (1 + z)(dV r /dr) 



(2) 



where A\q is the spontaneous emission coefficient for the 
transition (2.85 x 10~ 15 s -1 ), nm is the neutral hydro- 
gen density and dV r / dr is the gradient of the total radial 
velocity along the line of sight (with V r = V • n); on av- 
erage dV r /dr — H(z)/(1 + z). The neutral density can 
be expressed as 7Jhi = ^^(1 + 5b) (assuming the only 
baryon clement is hydrogen), when is the mean num- 
ber density of cosmic baryons, with a spatially varying 
overdensity 5b and xh is the fraction of neutral hydrogen 
(xh = 1 — x e where x e is the fraction of free electrons) . 
The 21-cm temperature is then (lijj : 



T 



T 21 (n,u) » T c 1- =M x 



"1 



(l + S b )(l + 5 XH )(l-\±^^\, (3) 



H(z) dr, 

where v is the peculiar velocity along the line of sight and 

h 



23 x H 



0.7 



0.02 



0.15 



10 



1/2 



mK. 



(4) 



The spin temperature is coupled to the hydrogen gas tem- 
perature (Tk) through the spin-flip transition, which can 



be excited by collisions or by the absorption of Lya pho- 
tons (Wouthuysen-Field effect; 0]) and we can write: 



1 



T 
22. 



Vtot 



I 



Vtot 



1 



T 
21 

T) 



K 



(5) 



where ytot = Ua + Vc is the sum of the radiative and 
collisional coupling coefficients. When the coupling to 
the gas temperature is negligible (e.g. y to t ~ 0), Ts ~ 
T 7 and there is no signal. On the other hand, for large 
Vtot, T s simply follows T K . The coupling coefficients are 
Va = ■lP a T±/27A w T 1 and y c = ^Ki- Q (T k ) n H T i ,/'iA w T 1 , 
where P a is the Lya scattering rate which is proportional 
to the Lya intensity, and ki_o is tabulated as a function 
of T fc Hill. 

To first order in the spin temperature perturbations 
related to collisions and radiation coupling, we then have 



1-S 



Vtot 



T s 1 

Va 



(1 + vtot)ytot 



ytot 
5 V „ 



T 



T~ 



Tk — T~, 



(1 + Vtot)Vtot Vc 



(6) 



The above perturbations contribute with different 
weights to the 21-cm temperature fluctuations depending 
on the redshift and we can roughly consider the following 
regimes of interest: 



1. z > 200 — Here, T s 
any signal. 



T~ and we do not observe 



2. 30 < z < 200 — The spin temperature approaches 
the gas temperature through collisions (y c ^> y a ). 
The gas is cooling adiabatically with Tk < T 1 and 
the signal is observed in absorption (see, Figure 1). 
Here, xh — 1 (no reionization). The 21-cm fluctu- 
ations are essentially sourced by perturbations in 
the gas density, temperature, perturbations in the 
collisional parameter and the radial velocity gra- 
dient [1, . These can be calculated from linear 
theory arguments without significant uncertainties 
in the astrophysics. Since the observations must 
be at very low radio frequencies, measurements are 
challenging and none of the planned interferome- 
ters focus on this regime yet. 

3. 2 hcat < z < ^sources — Initially, the coupling of 
Ts to T 7 dominates over the collisional coupling 
to Tk since the gas is rarefied through the expan- 
sion of the universe. As soon as the first galaxies 
begin to appear at ^sources, the Lya photons pro- 
duced by these sources couple Ts to Tk through 
the Wouthuysen-Field effect 0]. We then have 
Ts ~ Tk < T.y, xh — 1 at the beginning but could 
be changing as reionization begins, and y a ^ y c 
(collisions are only important at redshifts greater 
than 30 or so). Measurements of the brightness 
temperature in this region can provide information 
about the sources of Lya photons, the reionization 
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process, and the onset of heating, as well as cos- 
mology. We will focus on this regime since uncer- 
tainties remain as to when exactly Zhoat happened 
and the physics behind it. 



4. z re i ^ z z hcat - As the gravitational collapse con- 
tinues, and sources continue to reionize the Uni- 
verse, the gas temperature is heated above the 
CMB through a background of X-ray photons [l(| • 
The 21-cm temperature is then observed in emis- 
sion and Ts ~ Tr 3> T 7 . In this regime, pertur- 
bations are mostly dominated by variations to the 
Hydrogen neutral fraction, gas density and radial 
velocity gradient. This regime was already stud- 
ied in several papers 0, fla |: the projected 21-cm 
experiments will be able to follow the reionization 
history to a significant accuracy. If Ts T 1 be- 
fore patchincss of reionization is significant, then, 
one can use that epoch to test the cosmological 
model as the only relevant perturbation to 21-cm 
background is due to baryons (e.g. 6 X „ <C 6b). 
This is the case studied in Refs. [12|, [jjl, but it 
is unclear to what extent the Universe will remain 
mostly neutral (so as to produce a measurable 21- 
cm signal) while neither patchy nor cold (T s ~ T 7 ). 
It could also be that the redshift interval between 
Zhoat and z ro ; is small such that reionization rapidly 
follows heating; this will complicate the analysis as 
21-cm anisotropy measurement and foreground re- 
moval is then subject to prior assumptions on the 
astrophysics. 



Since there is a large uncertainty in the existence of 
a period where conditions allow simple theoretical and 
analytical calculations of 21-cm fluctuations we 
consider a more general case here. We study the situation 
where reionization began with a background of Lyman- 
a photons and gas heating followed subsequently well 
before the complete reionization of the Universe. This is 
equivalent to a combined case involving the third and the 
fourth regime from above. Note that our model ignores 
heating of gas within overdense regions well before the 
appearance of an X-ray background. While such heating 
is expected to be concentrated on small-scale structure as 
seen in simulations of Ref . [13] , this scenario could com- 
plicate both analytical predictions on the expected 21- 
cm signal as well as comparisons to observations. Fortu- 
nately, such complications are restricted to small angular 
scales and we focus mainly on large scales with multipoles 
below 9000. In the next Section, we present an analyti- 
cal formulation of 21-cm anisotropics in the presence of 
a Lyman-a field, when the universe is partially ionized 
and fluctuations are sourced not only by the density field, 
but also by variations to the ionizing fraction, velocities, 
among others. 



B. 21-cm background anisotropies: a general 



scenario 



In this paper, we will concentrate on the regime where 
Lya coupling cannot be ignored and Tk < T~ instead of 
Tk ^> T 7 . Combining equation ([3]) and equation ([6]), to 
first order in the perturbations, we can write: 



T 21 {h,v) « T 21 {v) [l + 6 b + 6 



+ 



1 + z dv s 



I + ytot y " H{z)dr 



(7) 



where we are already assuming that y a ~ ytot, as colli- 
sions are not important, and the spatially averaged 21-cm 
temperature is, 



T 21 {v)=T c 



ytot 



l 



ytot 



In. 

T K 



(8) 



Note that if the neutral fraction perturbations are large, 
one might need to consider also the term 6b6 x „, but since 
this is a second order term, we ignore it here. We can 
further relate the velocity perturbation to the cold dark 
matter density perturbation (6 C ), by using 



v(k, z) = i 



G 



k 2 (l + z)G 



6 c (k,z) 



(9) 



where v(k) is the Fourier mode of the peculiar velocity, 
G(z) is the dark matter perturbation growth function, 
G = dG/dt is the derivative with respect to conformal 
time, and we are assuming the case with no vorticity, as 
expected under the linear theory for perturbations. 
We finally have: 



T 2 i (ft, v) r>T 2 i(i/) 



6 XH (k,z) 



T, 



T K -T. 



d 3 k ^ 
(2tt) 3 

■6 TK (k,z) 



6 D (k)+6 b (k,z) 
1 



1 + 



n • k 



f(z)6 c (k,z) 



-6 Va (k, z) + 



(10) 



where 6 D is the Dirac delta function, f(z) is the sup- 
pression factor given by / = dlnG/dlna (a is the scale 
factor) and r = r(z)h, with r(z) = j Q cH^dz', the 
comoving distance to redshift z — v^ijv — 1. At the 
redshifts we are considering, it is safe to assume that 
6T K (k,z) = g(z)6b(k, z), where g + 1 is called the adia- 
batic index. As the strength of the coupling of the gas 
to the CMB decreases, the expansion of the gas becomes 
essentially adiabatic and we have g ~ 2/3. 

For an interferometer type experiment, the fundamen- 
tal quantity to be measured is the visibility which is di- 
rectly related to the spherical harmonic moment of the 
21-cm fluctuations @, Ojj], given by 



M= / dnY^T^n,^). 



(11) 
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The measured brightness temperature (T s ) corresponds 
to a convolution of the intrinsic brightness with some 
response function W v that characterizes the frequency 
resolution of the experiment: 



T s (fi^ ) = J dvW V0 (v)T 21 (n,v). 



(12) 



Using equation (fT0|) we can reexpress equation (fTTj) as 



ofm(fo) = (47r)i ; T 21 (^ ) 



(Sb(k, z ) 



+ S XH (k,z ) + 



T, 



■S TK (k,z) 



1 



1 + ytot 



S ya (k,z ))lp(k) - f(z )S c (k,z )Jp(k) 



(13) 



where we have dropped the zero mode, which only con- 
tributes to the monopole or the mean of the 21-cm bright- 
ness temperature; 21-cm interferometers are not sensitive 
to the mean temperature. The measurement is challeng- 
ing, but useful as the mean brightness temperature cap- 
tures the reionization history [22j. 
In equation (fT3|) . 



dv 
dr 
dv 
dr 



drW„ Q {v)j l (kr). 



(14) 



The spherical Bessel function is given by ji(kr) and 
j l (kr) is the second derivative of ji(kr) with respect to 
its argument. In deriving this form for the af m , we made 
use of the Rayleigh expansion for the plane wave given 
by 



and 



(h-kYe 



2„ikr 



l__d_ 

k 2 dr 2 



(15) 



(16) 



In order to simplify the calculation, we also assumed that 
the time dependent perturbations could be considered 
constant across the width of the window function so that 
they could be taken out of the integral over the Bessel 
functions. For large window functions, this simplification 
can be easily removed if needed. We can then write the 
angular power spectrum of T s (n, ^o), 



(17) 



in terms of the 3D power spectrum of the perturbed quan- 
tities in equation (|13| . Given two perturbations, a and 
d, the corresponding 3D power spectrum (P a d) is defined 
through: 



(5 a (k, Vl )5 d {k' ,V 2 )) 



(27T) 3 6 D (k + k')P ad (k^ 1 ,V 2 ). 

(18) 



To proceed, we need a model for the spatial inhomo- 
geneities in the Lya flux as well as the ionization frac- 
tion. The Lya flux inhomogeneities are associated with 
the locations of first ionizing sources. These sources 
are likely to be highly biased with respect to the den- 
sity field as densest regions are expected to form first 
stars. Analytical modeling in Ref. [2(| suggests the power 
spectrum of the Lya inhomogeneities to be of the form 
PyaVc, W = W 2 (k)P cc (k) where the window function can 
be considered as a damping term on the CDM fluctuation 
power spectrum P cc (fc, z). Thus, we take a parameterized 
model of the form: 



XKz) = b z a {k,z)e- k H ^P cc (k,z) 



(19) 



where b a (k, z) is a scale-dependent bias and Rl v corre- 
sponds to a characteristic size scale that corresponds to 
the Lya emitting patch given when the sources turned 
on (R Ly ~ lOOMpc). 

For the ionization fraction field, we follow Ref. (l6l.[28l| 
and take a model following: 



Px H x H (^i z ) 



bl H (k,z)e- k2R *»P c , 



(k,z) . 



(20) 



where now the characteristic size scale R XH is that of 
reionization patches. Analytical models suggests that 
this scale is of order a few Mpc and grows to about 100 
Mpc towards the end of reionization as bubbles begin to 
overlap. Here, we assume 50 Mpc and 6 Mpc as two pre- 
ferred size scales for reionization bubble sizes correspond- 
ing to redshift bins of 7.5 to 8.5 and 8.5 to 9.5, respec- 
tively. We also consider an alternative scenario where 
bubble sizes may be lower than suggested by current 
models with values of 1.0 Mpc and 0.5 Mpc in redshift 
bins of 7.5 to 8.5 and 8.5 to 9.5, respectively. The bias 
term b XH corresponds to (xh — l)/xjjb where b is the bias 
used in Ref. [l6l |. In our calculations, we assume that the 
baryon density field is a biased tracer of the dark matter 
density field and write Pbb(k,z) = b 2 (k, z)P cc (k, z). To 
simplify our analytical calculation, we assume that the 
cross-correlation spectra between each of the fields con- 
sidered can be described with perfect correlations, such 
that Pad = ±\/ Paa.Pd.d- Note that 8 XH is ant icor related 
with S c . Any departures from perfect correlations will 
only lead to additional parameters and a further degra- 
dation of parameter estimates than suggested here. 

Putting all terms together, we finally have: 



Cf(v 1 ,v 2 )=T2i(v 1 )T2i(y 2 )- / k 2 dk x 



F(k,u 1 )F(k ) u 2 )iP(k)I^(k) 



+f(z 1 )f(z 2 )J?(k)Jp>(k) - F(k,^)f(z 2 )ir i (k)J^(k) 



-F{k,v 2 )f{ Zl )I^{k)J^{k) 



(21) 
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where 
F(k,v\z) 



fa(k,z) + b XH (k,z)e 
T 

s — :L 77r9{z)h{k 1 z) 
±K — i-i 



H' 2 {k,z)+ (22) 
L/2 



b a (k, z) 



1 + Vto 



Note that, if the visibilities follow a Gaussian distribu- 
tion, and as the radial width becomes smaller, the above 
angular cross-frequency power spectrum, contains all the 
three-dimensional information necessary to describe the 
21-cm temperature perturbations. This calculation as- 
sumes full cosmic evolution so that we don't need to 
worry about extra corrections coming from the pertur- 
bations time dependence as described in Ref . [2!| . More- 
over, the Alcock-Paczynski effect |30j is implicitly in- 
cluded in the calculation [3l[, though the line of sight 
radial modes are smoothed over a bin width of 0.1 MHz 
so we do not necessarily have information from radial 
modes at scales smaller than the three-dimensional power 
spectrum here. In the presence of foregrounds, and with 
interferometers that are mostly sensitive to the modes 
along the line of sight alone, the task of separating the 
three-dimensional power spectrum to various contribu- 
tions depen ding on the line of sight angle [HJ is extremely 
challenging [13| . 



III. FIDUCIAL MODEL AND PARAMETERS 

A. Pre-reionization: In the presence of Lyman-a 
photons 

The main objective of this paper is to determine how 
well one will be able to measure the cosmological param- 
eters and the astrophysical parameters related with the 
pre-reionization process when b XH — 0, Tk < T 7 and 
Lya coupling of T5 to Tk is dominating. Following the 
description in section Hi A[ we will assume that this re- 
gion corresponds to the redshift range between z = 15 
and z = 25. This is a reasonable assumption, specially 
on the high end, since collisional coupling should only be- 
come important for z > 40, though the exact end points 
of this region do not affect the main conclusions of our 
analysis. Our model Jgnores complex situations such as 
the one seen in Ref. [271 ] . where the 21-cm signal appears 
in emission in overdense regions well before IGM heat- 
ing by an X-ray background. However, this emission is 
restricted to small halos corresponding to structures at 
multipoles larger than 4 x 10 4 for z ~ 17 and next gener- 
ation experiments are not expected to have the required 
sensitivity to measure the power spectrum at such small 
angular scales. At larger angular scales the fluctuations 
are determined primarily by the absorption and the emis- 
sion signal is unlikely to impact our analysis. Since the 
region we "observe" corresponds to the pre-reionization 
scenario, we set b XH = and ignore the presence of par- 
tially ionized bubbles. We could, however, consider the 
most general scenario in equation I22[ including the ion- 
ization fraction term, the collisional coupling and adding 
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FIG. 1: Temperature evolution: CMB (solid line), gas 
(dashed) and Spin (dotted). The gas is assumed to cool adi- 
abatically for z < 145 so that Tk ex (1 + z ) 2 ■ Note that for 
z > 30 collisions will become important, so that Ts ~ Tk 
again. 



a few more parameters to account for the heating of the 
gas above the CMB. This would however, make it ex- 
tremely hard to distinguish between the different param- 
eters just by using the anisotropy power spectrum alone, 
function of the frequency. 

Note however that there should be a well defined sig- 
nature for this epoch: the average 21-cm temperature is 
smaller than the CMB temperature, T21 < in equation 
[7] (e.g. the signal is observed in absorption). For higher 
redshifts T$ will increase towards T 7 before collisional 
coupling dominates, so that T 2 i ~ 0. At lower redshifts, 
the gas temperature will increase above the CMB tem- 
perature so that T21 > 0. While T21 < 0, it is safe 
to neglect the perturbations in the ionization fraction so 
that we can simplify the calculation, though it is more 
complex than some of the discussions in the literature so 
far. 

At these redshifts, the gas is cooling adiabatically with 
T K (x (1 + zf and g(z) ~ 2/3 in equation [H The 
gas temperature should equal the CMB temperature at 
z ~ 145 when T 7 « 397.850X. For y a we follow Ref. [H 
(see figure 5 in their paper; also, 3). Figure [1] shows the 
evolution of T 7 , Tk and T$ in our model. 

As stated above, eventually the gas will be heated 
above the CMB so that the signal is observed in emission 
and reionization starts, but we assume this won't happen 
until z < 15. Even if there is some heating of the gas at 
the low end of the redshift interval, the parameterization 
we will be using can still account for that as long as the 
bias associated with the gas temperature perturbations 
is independent of scale. The "astrophysical parameters" 
are Rl v and because of parameter degeneracies, we com- 
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bine the other unknowns in equation (|2ip into: 



70) 
(3{k,z) 
a(k 7 z) 



Vtot 



1 + Vtot 

bbik, z) 
b a (k,z) 



1 - 



Ik 

T 



K 

T-. 



T K - T- 



■9{z) 



1 



Vtot 



(23) 
(24) 
(25) 



Note that, although we are neglecting perturbations in 
the neutral fraction b XH = for the redshift range con- 
sidered, we can still easily accommodate for a smooth 
reionization component by allowing xh to be a free pa- 
rameter (just in case reionization starts when Tk < T 7 ). 
This parameter will however, be degenerate with 7(2) 
through T c (see equation [SJ • 

We also consider the following cosmological parame- 
ters: il m h 2 , Ob/i 2 , f^A, n s and assume a fiat-cosmological 
model. For the primordial power spectrum, we use: 



p 



rn s +3 



(26) 



where the power spectrum normalization, Sh is a free 
parameter. This is also degenerate with xh and 7(2) as 
they all enter in as a product in equation (|2"Tj) . When 
combined with CMB, this degeneracy can be broken and 
T c dependence is determined through other cosmologi- 
cal parameters allowing the product x~h"/(z) to be de- 
termined. It is unlikely that one can independently de- 
termine the neutral/ionizing fraction and the coupling 
strength of Lya background, j/tot ~ Da from 21-cm 
anisotropy observations (unless it is safe to assume that 
Ts >> T 7 ). This will complicate the attempts to deter- 
mine the reionization history of the Universe from 21-cm 
observations alone. 

Table U shows the fiducial values assumed for the time 
independent parameters in our model. The cosmolog- 
ical parameters correspond to the best fit values from 
the WMAP 3-year data [23j. Note that 6h is chosen so 
that the fiducial model gives o~% — 0.74. Figures [21a) 
and (b) show the expected evolution of the time depen- 
dent parameters for the redshift range we consider in the 
analysis. 

Figure [3] shows the 21-cm angular power spectrum, cal- 
culated using the above model, as a function of I, for the 
largest and smallest frequency we consider in our range 
(55 MHz and 85 MHz). As the frequency increases (red- 
shift decreases), the power also increases. Also shown is 
the separate contribution to the total power from the ve- 
locity perturbations, the density perturbations and the 
velocity-density correlations. Note that the power spec- 
trum is just the sum of these three components. Since 
the velocity perturbations are a direct probe of the dark 
matter perturbations, if we could measure this compo- 
nent it would be possible to make a clean measurement 
of the cosmological parameters without the "contamina- 
tion" from the other astrophysical parameters 20] . Fig- 
ure 2] shows the same quantities but as a function of fre- 
quency. In order to parameterize the redshift evolution of 
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FIG. 2: Top: Evolution of 7(2). Note that 7 is always neg- 
ative since Tk < T 7 and the moduli of 7 is increasing with 
decreasing redshift due to the cooling of the gas, which gives 
a boost to the 21-cm signal. Bottom: Evolution of f3{z) (solid 
line) and ct(z) (dashed). 



the 7, j3, a parameters, we considered small redshift bins 
of Az ~ 2 and assumed these parameters were constant 
in that bin. 



B. During reionization: when Tk » T- t . 

In order to compare our results with prior studies fl2l . 
[l3l ] , we also considered in our analysis the case when Ts = 
T K » Try, that should be valid at redshifts below z = 15 
in our model. Note however that, although by observing 
the signal in emission we know that Tk > T 7 , there is no 
guarantee that we can indeed ignore the perturbations in 
the gas temperature and in the Lya coupling. If Tk >> 
T 7 , then the calculation is significantly simplified and 
21-cm fluctuations are primarily sourced by the density 
perturbations in the gas field and the ionization fraction. 

As reionization proceeds, while it will be safer to as- 
sume Ts » Tj, one cannot ignore the perturbations 
in the ionization fraction which might dominate the sig- 



8 



TABLE I: Fiducial parameters used in our analysis. 





^ 77? h 


Q b h 2 


Q A n s 8 H x 10 5 


Rhy 


Fiducial values 


0.127 


0.0223 


0.762 0.951 6.229 


100 Mpc 



The cosmological parameters based on the WMAP 3- year data |2.'J . We take a fixed size for the inhomogeneities in the Lyman-a 
radiation density field, with a size scale of Rhy (see, equation ll9ll ; this is only an approximation and the true scenario with evolving 
inhomogeneities can only complicate the analysis than suggested here. In addition to Rl v , astrophysical parameters include a(z), j3(z), 
~f(z), and Xfj(z) and we set their values as a function of the redshift bin of interest (see, Figure 2, for example). 




FIG. 3: The 21-cm power spectrum (solid line) and the den- 
sity (dashed), velocity (dotted) and density- velocity (dot- 
dashed) contributions. Note that the total power is just the 
sum of these three contributions. Bottom (blue) lines corre- 
spond to v — 55 MHz (z ~ 24.8) and top (red) lines corre- 
spond to v = 85 MHz (z ~ 15.7). 
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estimating, for example, cosmological information in the 
21-cm background [131 ]. There are large uncertainties as 
to when the gas is heated to significantly higher temper- 
atures above that of the CMB, since one does not have a 
good understanding of the primordial X-ray background, 
or the sources that contribute to this background, that 
are responsible for the heating. It could be that heat- 
ing followed the onset of full reionization, in which case 
the 21-cm signal will be small as the fraction of neutral 
Hydrogen is not significant. We take this idealized case 
simply as an example since we can compare our results to 
previous estimates in the literature. We do not, however, 
neither motivate this simple description nor suggest cos- 
mological precisions within this description as the true 
scenario. To be safe, one must include all sources of per- 
turbations when model fitting the 21-cm background. 



IV. THE MEASUREMENTS 

In this section, we will briefly describe our assumptions 
of instrumental noise associated with upcoming interfer- 
ometers as well as the foreground signals at low radio 
frequencies that must be cleaned to extract the cosmolog- 
ical 21-cm signatures. In addition to several experimen- 
tal setups, we will also outline a more general reference 
experiment that is optimized for anisotropy observations. 



Noise 



The noise in the angular power spectrum is given by 



c, 



N 



A 2 T< 



St) Si 



A 



dish 



(Au) 2 
Avt v (l)' 



(27) 



FIG. 4: The 21-cm power spectrum (solid line) and the den- 
sity (dashed), velocity (dotted) and density- velocity (dot- 
dashed) contributions as a function of frequency (in MHz). 
Bottom lines (thin) correspond to I — 10 and top lines (thick) 
correspond to I = 6700. 



nal at lower redshifts @, [l6j]. The latter source of in- 
homogeneity is generally ignored in the literature when 



where Au is the width of the primary beam and is set by 
the size of the dishes d, Au ~ d/ A. The observation time 
per visibility, t v , usually depends on the the multipole I 
being observed and can be approximately expressed as: 



t V 



t N(l)Al 
2^1 : 



(28) 



where to is the total observation time and N(l) corre- 
sponds to the number of baselines contributing to the 
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multipole I. For uniform Fourier coverage, the number of 
baselines must grow linearly with I and we have 



N(l) 



2wlAl 
tt! 2 

max 



N h 



(29) 



Note that Al = 2nAu, JV& is the total number of baselines 
and Imax corresponds to the largest visibility observed 
within this uniform coverage. We finally have 



2 12 ' 

cover 1 max 



(30) 



where f cover is the fraction of the total area covered by 
the dishes: 



fc 



_ NdishAdish Ndish(Aiy 



A 



total 



I 2 

rn a x 



(31) 



For realistic arrays, the number of antennas does not in- 
crease linearly with the size of the baselines, so that the 
Fourier coverage is not uniform (e.g. N(l) does not grow 
linearly with I) and the noise power spectrum will vary 
with I. Typically, the array will have a core at the cen- 
ter where the antennas are closely packed (large covering 
fraction) and a more dilute configuration for larger base- 
lines (smaller covering fraction at higher Us). We will 
take into consideration the possible effects of the geom- 
etry of the array by writing the noise power spectrum 
as 



T 2 



(27T) 



Aut i 2 f 2 (iy 



where, for uniform coverage, 

hi) = fco 



i 



(32) 



(33) 



The exact antenna distribution has not been decided 
for any of the interferometers we consider here. To pro- 
ceed, we will assume that the Fourier coverage is uniform 
within the core of the array and for each of the possible 
outer annuli, but allow for different coverage densities 
for each of these regions. If, for instance, the density of 
antennas goes as ~ r~ 2 within the core, then the actual 
noise will be lower than what is assumed here close to the 
center and larger away from the center. This small differ- 
ence should cause very little impact on the final results. 
We will consider several experimental specifications here 
related to upcoming low-frequency interferometers that 
are either under construction or planned. 

First among the upcoming interferometers is the Low 
Frequency Array (LOFAR) that will have a total collect- 
ing area of about 2 x 10 5 m 2 with approximately 25% 
of that area concentrated in a compact core of radius 
R = 1 km and 50% of the antennae within a 6 km ra- 
dius. For the core we have /(/) ~ 0.016l max /l and for 
the outer region f(l) ~ 0.001l max /l. Note that l max is 
given by 



4-7T- 



R 



(34) 



which corresponds to l m ax ~ 5500(f/130MHz) for the 
inner core of LOFAR. 

The MWA array consists of 500 antennas distributed 
within a 0.75 km radius [33j. The total collecting area 
is approximately 8000m 2 so that f(l) ~ Q.0045l ma x/l- 
There is also a proposed experiment based on an ex- 
panded MWA configuration (MWA5000) with a total of 
5000 antennas of the same type. The antennas are ex- 
pected to be distributed over an area of 1.5 km radius, 
with approximately 80% of the antennas within a radius 
of 0.75 km, so that /(/) ~ 0.036l max /l for the inner ra- 
dius and f(l) ~ 0.011l max /l for the outer core. 

A more ambitious project is the SKA that will have 
an effective area of ~ 10 6 m 2 . Since details of the array 
are still not finalized, we will assume that (20%, 50%, 
55%) of the baselines will be within a (0.5, 3, 6) km 
radius. The correlation length in the visibility space, 
Al = 2wd/X, is set by the diameter of the "dishes", d, 
which is approximately 4 m for MWA and MWA5000, 35 
m for LOFAR and 10 m for SKA. Note that the diameter 
of the antennas does not necessarily set the lower limit on 
the baseline length since the antennas might not be closed 
packed together. This is especially true for LOFAR where 
the minimum baseline is 100 m giving l m in ~ 270 at 
z = 130MHz. Table ITTl summarizes the specifications for 
each of the experiments considered. 



B. Reference experiment 

In order to probe the 21-cm signal on the relevant 
scales (I ~ 1000) at the high redshifts we are considering, 
we need baselines of the order of 8 km. This is outside the 
scope of most of the experiments described above, except 
SKA. The sensitivity to the signal is also complicated by 
the fact that the system temperature is completely dom- 
inated by the huge sky temperature. In our calculations 
we used T sky oc (1 + z) 2 6 (with T sky ~ 1500^ at z = 15). 
As discussed above, for SKA we assumed the minimum 
baseline length to be d=10m which will set a beam size 
of Al ~ 18 at z=15 and Al ~ 11 at z=25. 

Taking into account the above variations, we assumed 
it was possible to generate an effective beam with fixed 
size Al ~ 20. The constancy of this beam will be es- 
sential to properly remove the foregrounds. Ignoring for 
the moment the residuals due to foregrounds, we show in 
figure [5] the expected errors in the measurement of the 
power spectrum when using an experiment like SKA. We 
assume a total of 2000 hours of observation for 2 regions 
of the sky and a frequency resolution of Av ~ 0.1 MHz. 
Note that the field of view of these experiments is fully 
taken into account in the error calculation (e.g. sample 
variance is included). 

We see that the expected error for SKA is well above 
the signal at z ~ 25. This is essentially due to the large 
noise temperature associated with the sky at these low 
radio frequencies. For comparison, we also show the 
error for a fiducial experiment that has 100% coverage 
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TABLE II: Parameters adopted for the experiments. 











Radius 


LUV CI CLtiC i.1 dL LIU 11 


Array 


effective 


area 


antenna diameter Rl, R2, R3... (km) 


fl, f2, £3... 


LOFAR 


2 x 10 5 


m 2 


35 m 


1.0, 6.0 


0.016, 0.001 


MWA 


8 x 10 3 


m 2 


4 m 


0.75 


0.0045 


MWA5000 


8 x 10 4 


m 2 


4 m 


0.75, 1.5 


0.036, 0.011 


SKA 


1 x 10 6 


m 2 


10 m 


0.5, 3.0, 6.0 


0.25, 0.018, 0.0049 


SKAb 


1 x 10 6 


m 2 


10 m 


0.5, 4.0 


0.38, 0.018 




l 

FIG. 5: The 21-cm power spectrum (solid lines) and the ex- 
pected errors. Thin red lines show the results at z ~ 15 while 
thick blue lines show for z ~ 25. Dashed lines shows the ex- 
pected errors for SKA, dotted lines for a fiducial experiment 
with 100% coverage and dot-dashed lines for another SKA 
type experiment (our reference experiment - SKAb). 

fraction up to a radius of ~ 4 km. Even for this case, 
the only way to measure the power spectrum is to in- 
crease the observation time (or the number of observed 
sky patches). Note however that these results are for 
only one frequency channel with 0.1 MHz resolution. In 
principle we could pack together many more channels in 
our analysis as long as cosmic evolution can be neglected. 
Moreover, there will also be sonic extra information from 
the cross-frequency correlations. We therefore decided 
to proceed our analysis using an SKA type experiment 
with a few changes to improve the signal to noise at very 
low frequencies. We assumed the same collecting area as 
SKA, but with 30% of it within a 0.5 km radius and 90% 
of the antennas within a 4 km radius. We call this array 
configuration SKAb. We used a total of 3000 hours of 
observation time for 2 locations in the sky. 

C. Foregrounds 

We take into account the expected foregrounds at these 
frequencies, following Ref. [lj|- Table HLT1 shows the fidu- 



cial values for the free foreground parameters used in our 
analysis. We model the foregrounds as power laws in 
both I and v so that: 

Cf{v, v) = A t {l^/lf'{u f /u) 2& \ (35) 

where Vf = 130MHz. We also expect the foregrounds to 
be highly coherent, i.e. I™(^i, v?) ~ 1 where 

If K u 2 ) = CfK ^)/^CfK^)CfKzA0. (36) 

Although we consider other parameterizations of our 
starting point is 

Zf K v % ) = exp [- \oi( Vl / V2 )/2/g] , (37) 

which for the frequency range and the values of consid- 
ered, can be written as J\ % {yi, v-i) ~ 1— log 2 (i/i/^2)/2/£?i- 
This form captures departures from perfect correlations. 
This decorrelation will also lead to a departure from a 
power-law in the frequency dependence of the foreground 
brightness temperature, with a varying as a function of 
the sky position. We refer the reader to Ref. [IB] for more 
details. 

In addition to foregrounds, other sources of contam- 
ination at low-frequency 21-cm observations include in- 
terferences and systematics. We ignore such subtleties 
here. Ignored here is also the confusion that might arise 
from radio recombination lines [34[, which do not pro- 
duce smooth foregrounds in the frequency space. Radio 
recombination lines may require the removal of certain 
frequency channels, provided that the presence of radio 
recombination lines can be identified at high frequencies 
above those related to 21-cm observations along a given 
line of sight. 



V. RESULTS ON PARAMETER 
MEASUREMENTS 

In order to establish the extent to which both astro- 
physical and cosmological parameters can be extracted 
from the data, we make use of the Fisher matrix formal- 
ism. We work in the angular multipole space with power 
spectra binned in the frequency resolution of Ais = 0.1 
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TABLE III: Fiducial parameters for foregrounds. 





A(mK 2 ) 


P 


Q 


e 


extragalactic point sources 


10.0 


1.1 


2.07 


1.0 


extragalactic free-free 


0.014 


1.0 


2.10 


35 


galactic synchrotron 


700 


2.4 


2.80 


4.0 


galactic free-free 


0.088 


3.0 


2.15 


35 



The amplitude of the extragalactic point sources is calculated 
assuming a flux cut (S c ) of 3 jiJy which is typical of SKA. This 
amplitude roughly scales as A ~ Sj?' 5 . 



MHz throughout our calculations. Although the experi- 
ments can go well below this value, there is little infor- 
mation to be gained from using better resolution due to 
the correlations of the signal between different adjacent 
frequency bins |16j |. 

Instead of the three-dimensional power spectrum of the 
brightness temperature, Pr(k), we make use of the angu- 
lar power spectra as the observable. The measurement of 
Pr(k) in Ref. [l3| is motivated by the fact that one can 
separate it into terms with different multiplicative factors 
associated with the line-of-sight angle [321] . We do not 
implement such a possibility here as the upcoming inter- 
ferometers arc expected to provide mostly information on 
modes along the line of sight with no significant informa- 
tion on the modes perpendicular to the line of sight [l3T ] , 
which are needed to separate the terms related to, say, 
density perturbations from the velocity term. Our for- 
mulation includes the information from velocity pertur- 
bations and other sources of directional-dependent inho- 
mogeneities. While we do not consider the reconstruction 
of three-dimensional P(k, z) from two-dimensional Ci an- 
gular power spectra, by binning the measurements over 
a bin of 0.1 MHz and by including the cross-correlation 
information between different frequency bins, C;(i^,za,-) 
as the measurement, we capture information on the ra- 
dial modes of three-dimensional clustering, but averaged 
over a comoving distance of L ~ 1.7-*/ (1 + z)/lQ that 
corresponds to the used frequency bin of 0.1 MHz; Thus, 
the only difference between this analysis and some of the 
suggestions in the literature is that we do not have in- 
formation on radial modes below this length scale. Our 
smoothing scale, however, is adequate enough to capture 
information from, for example, the velocity fields and 
the correlation of velocities across adjacent radial bins. 
We believe this is an adequate and a conservative ap- 
proach to take when analyzing 21-cm information, given 
the lack of detailed information related to foregrounds at 
low radio frequencies and how they may impact detailed 
three-dimensional power spectrum measurement. 

The covariance at a given multipole between two fre- 
quency bins is 

Cl = C{(yi, Vj ) + CUvuVj) + Cf , (38) 



from which we write the necessary Fisher matrix as 

when p a and p are two parameters of interest. In here, 
N(l) is the number of independent modes measured by 
an experiment. The inverse of the Fisher matrix pro- 
vides the minimum error to which parameters can be 
established from the data. 

We calculate the Fisher matrix and then the expected 
errors in the cosmological and astrophysical parameters 
assuming we have data (visibilities) from v = 55MHz 
to v = 85MHz (z = 15.7 and z = 24.8). We also take 
into account the expected foregrounds at these frequen- 
cies and use a variety of experimental descriptions, with 
parameters in Table II. Basically, foregrounds can be re- 
moved due to their smoothness (high correlation) across 
the whole frequency range, while the signal is oscillating 
widely from one bin to another; as the correlation of the 
signal across frequency bins increases it becomes harder 
to identify and remove the foregrounds. 

Figure [6] shows the signal frequency cross-correlations 
for a few multipolcs. We see that the signal is highly 
correlated for Av = 0.1 MHz, unless we consider scales 
of I > 10000. Moreover, for large scales (small I), we 
need to consider a large frequency range in order to dis- 
tinguish the foregrounds from the signal. Therefore we 
also need to be careful with the frequency range one con- 
siders in the analysis. From a practical point of view it 
will be easier to do the analysis for small intervals (say 8 
MHz) between the 55 MHz and 85 MHz and then com- 
bine the results. However, if the foregrounds are not 
properly removed in these small intervals, the final result 
will be worse than the residual foreground level with the 
full range. These issues require a more detailed analysis 
of foregrounds, which will only happen after first obser- 
vations of the low frequency radio sky with instruments 
such as MWA and LOFAR. 

In order to see the effect of the interval size, table 
IIVI shows the combined errors in the parameters from 
a Fisher analysis using a 5 MHz, 10MHz and 20 MHz 
interval around v = 70 MHz (z « 19.3). We considered 
a, (3 and 7 to be constant across these intervals, as well 
as the instrumental noise. Note also that we made the 
simplification that the foreground parameters are known 
a priori. If foregrounds were completely removed for any 
of the intervals, then the errors in the parameters, er, 
would just scale as cr/y^n), were n = 2(4) for the 10(5) 
MHz interval. This is clearly not the case, even for the 
10 MHz interval, so we will opt to use the full frequency 
range in our analysis. 

A. When T s = T K < T 7 and with Ly - a Coupling 

We did a Fisher matrix analysis taking into account 
the full three dimensional "angular" power spectrum 
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TABLE IV: Errors on cosmological and "astrophysical" parameters for three frequency intervals. 
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a 


R Ly (Mpc) 




Qf,h 2 




n s 


NS 


Values 


-3.13 


0.223 


0.48 


100 


0.127 


0.0223 


0.76 


0.951 




cr (20 MHz) 


26.48 


0.036 


2.11 


281 


0.773 


0.316 


0.10 


0.299 


4.6 


cr (10 MHz) 


210.87 


0.060 


20.45 


1673 


6.045 


2.519 


0.19 


1.275 


25.5 


cr (5 MHz) 


666.52 


0.110 


507.78 


18039 


18.755 


7.902 


0.39 


3.414 


122.1 



NS = l/AfJ^i kTu7n ' ' s * ne va l ue °f the parameter i and n = 2(4) for the 10(5) MHz interval. If foregrounds were completely 
removed then errors should scale as a/yfn and NS should be roughly constant. 
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FIG. 6: The 21-cm correlations at v = 65 MHz (z ~ 20.9), 
for I = 5 (dashed line), I = 270 (dotted line) and I = 10000 
(solid line). 



Ci{v x ,v 2 ), from v = 55MHz to v = 85MHz with Av = 
0.1MHz (or from z = 24.8 to z = 15.7). Note that this 
high redshift regime is unlikely to be probed with the 
first-generation 21-cm interferometers and thus we con- 
sider observations with our fiducial experiment, SKAb, 
which has the best noise temperatures at low frequen- 
cies of all experiments considered here. We parame- 
terize a(z), (3(z) and 7(2:) to be constant within bins 
of Az — 2MHz. All other time dependent factors in 
Ci{v\,V2) (e.g. the growth of perturbations) are allowed 
to fully evolve within these redshifts bins. We also in- 
clude foreground information as described in the text and 
when estimating errors in cosmological and astrophysical 
parameters, we marginalize over the uncertainty in the 
Foreground model. Figures [3 [8] and [9] show the expected 
errors for these time dependent parameters using SKAb. 

In Table [Vj we also show the forecasted 1 — a un- 
certainties for the cosmological parameters using SKAb 
with and without information on cosmological parame- 
ters from Planck, which will break the degeneracy with 
the amplitude. 

Table PVTl shows the expected correlations between cos- 
mological parameters as well as the parameters that cap- 
ture astrophysics; note that we have taken xh = 1 here. 



As listed in Table IVTl there are significant correlations be- 
tween parameters 74 and cosmological parameters of n s , 
f2b/i 2 , and Q, m h 2 . While the combination with Planck 
allows these correlations to be broken so that parameters 
such as "fi are independently established (see the improve- 
ment in errors in Figures [7] and [5]) , there is no significant 
information on cosmological parameters themselves from 
21-cm observations. Note that this is for the more gen- 
eral case of anisotropies in the 21-cm background at high 
redshifts; for more limited scenarios, as described below, 
we do find additional improvements. Figure [10] shows a 
few error ellipses for some of the more correlated param- 
eters in table IVI1 before and after adding Planck priors 
(involving parameter combinations of jiXH(zi) vs. cos- 
mological parameters such as fi m /i 2 and fli,h 2 ). Again 
we see that with 21-cm data alone the astrophysical pa- 
rameters are degenerate with cosmological ones, but with 
CMB information from Planck, cosmological parameters 
are pinned down allowing a measurement of parameters 
related to astrophysics. 

While the improvement in cosmological parameter es- 
timates is not significant at these high redshifts, 21-cm 
observations with an instrument such as SKAb provides 
high signal-to-noise ratio estimates of astrophysical pa- 
rameters such as 7i and Pi as a function of the redshift 
bin i. These parameters capture the thermal state of the 
gas relative to that of the CMB. Establishing (3 < 1 alone 
would be sufficient to show that the we are in a regime 
where Tk % T^. 



B. When Tk » T~ t and no fluctuations in the 
ionization fraction. 

To consider the extent to which parameters can 
be measured with the first-generation experiments, we 
moved to slightly lower redshifts and assume the scenario 
when Ts >> T 7 . In this regime, fluctuations in 21-cm 
brightness temperature comes from the inhomogencitics 
in the gas density and the ionization fraction. We first 
consider the case that b XH ~ such that fluctuations in 
the neutral/ionization fraction are not important; the an- 
gular power spectrum of 21-cm anisotropies is illustrated 
in Fig. QT] Table IVIII shows the forecasted 1 — a uncer- 
tainties for several experiments when using data between 
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TABLE V: Forecasted 1-a uncertainties on cosmological and astrophysical parameters when Ts < T 





R Ly (Mpc) 






n s 


5 H x 10 5 


Values 


100 


0.127 


0.0223 0.76 


0.951 


6.229 


SKAb 


95 


0.275 


0.1127 0.02 


0.075 




Planck 




0.0023 


0.00017 0.011 


0.0047 


0.03 


SKAb + Planck 


39 


0.0020 


0.00017 0.0096 


0.0044 


0.03 



Note: we assumed a total of 3000 hours of observation on two places in the sky and used a frequency interval between 55 MHz and 85 
MHz (we allowed for full cosmological evolution within this interval and also considered the expected foregrounds at these frequencies 

with models from 



TABLE VI: Parameter correlations using SKAb / SKAb + Planck when Ts = Tk J$ T y and with Lyman-a photons present. 





7i 


72 


73 


74 


75 


pi 


Pi 


P3 


Pi 


Ps 


7i 


1.00 




















72 


1.00/0.80 


1.00 


















73 


1.00/0.24 


1.00/0.17 


1.00 
















74 


0.81/0.03 


0.81/0.03 


0.81/-0.07 


1.00 














75 


0.09/0.00 


0.09/0.00 


0.09/-0.01 


0.05/-0.04 


1.00 












pi 


0.03/0.51 


0.03/0.25 


0.03/0.08 


0.02/0.00 


0.00/0.00 


1.00 










Pi 


0.05/0.10 


0.06/0.64 


0.04/-0.03 


0.04/0.00 


0.00/0.00 


0.10/-0.02 


1.00 








P3 


0.03/0.02 


0.03/-0.02 


0.12/0.95 


-0.02/-0.08 


-0.01/-0.01 


0.02/0.01 


-0.05/-0.06 


1.00 






Pi 


0.01/0.00 


0.01/0.00 


0.00/-0.08 


0.57/0.97 


-0.03/-0.03 


0.00/-0.01 


-0.00/-0.01 


-0.09/-0.09 


1.00 




Ps 


-0.01/0.00 


-0.01/0.00 


-0.01/-0.01 


-0.04/-0.05 


0.97/0.98 


0.00/0.00 


0.00/0.00 


-0.01/-0.01 


-0.04/-0.04 


1.00 


Ql 


-0.67/0.02 


-0.67/0.02 


-0.67/0.02 


-0.53/0.03 


-0.07/-0.01 


0.04/0.03 


i~\ c\-\ 1 f\ r\ a 

0.01/0.04 


0.01/0.03 


0.02/0.03 


0.00/-0.01 


OL2 


-0.36/0.00 


-0.36/0.04 


-0.35/0.03 


-0.27/0.03 


-0.05/-0.02 


0.03/0.00 


0.06/0.08 


0.03/0.04 


0.03/0.04 


-0.01/-0.02 


a 3 


-0.04/0.00 


-0.04/-0.01 


-0.03/0.17 


0.00/0.05 


-0.02/-0.01 


0.03/-0.01 


0.00/0.00 


0.15/0.15 


0.04/0.04 


-0.02/-0.02 


Q 4 


0.04/0.00 


0.04/-0.01 


0.03/-0.05 


0.17/0.24 


0.02/0.02 


0.00/-0.01 


0.00/0.00 


-0.04/-0.04 


0.21/0.21 


0.01/0.01 


as 


0.01/0.00 


0.01/0.00 


0.01/-0.01 


-0.03/-0.07 


0.41/0.41 


0.00/0.00 


0.00/0.00 


-0.01/-0.01 


-0.05/-0.05 


0.38/0.38 


RLy 


-0.86/-0.02 


-0.86/0.00 


-0.85/0.01 


-0.69/0.03 


-0.09/-0.02 


0.04/-0.01 


-0.02/0.04 


-0.01/0.04 


0.01/0.04 


0.00/-0.02 


^ ?? i h 


0.99/-0.02 


0.99/-0.03 


0.99/-0.02 


0.81/0.00 


0.09/0.00 


-0.02/0.11 


0.04/0.02 


0.03/0.00 


0.01/0.00 


-0.01/0.00 


n b h 2 


1.00/0.08 


1.00/0.06 


0.99/0.02 


0.81/0.00 


0.09/0.00 


-0.01/-0.02 


0.04/0.00 


0.03/0.00 


0.01/0.00 


-0.01/0.00 


n A 


-0.02/-0.02 


-0.02/-0.03 


-0.02/-0.02 


-0.02/-0.03 


-0.01/0.00 


-0.59/-0.34 


-0.23/-0.11 


-0.04/-0.02 


-0.01/-0.01 


0.00/0.00 


n 3 


-0.79/0.23 


-0.79/0.18 


-0.79/0.05 


-0.64/0.01 


-0.08/0.00 


0.16/0.11 


-0.01/0.02 


-0.01/0.00 


-0.01/0.00 


0.00/0.00 


A 


-/0.07 


-/0.06 


-/0.02 


-/0.00 


-/0.00 


-/0.00 


-/0.00 


-/0.00 


-/0.00 


-/0.00 




CYl 






a 4 




Rhy 


in h 




Qa 


n s A 


ai 


1.00 




















a.2 


0.78/0.78 


1.00 


















a-i 


0.38/0.47 


0.56/0.56 


1.00 
















Ct4 


0.06/0.11 


0.16/0.17 


0.44/0.44 


1.00 














«5 


0.05/0.07 


0.10/0.11 


0.28/0.28 


0.39/0.39 


1.00 












Rl v 


0.85/0.85 


0.71/0.90 


0.34/0.60 


0.06/0.17 


0.03/0.08 


1.00 










^ m h 


-0.67/-0.01 


-0.38/-0.08 


-0.07/-0.08 


0.03/-0.03 


0.00/-0.01 


-0.88/-0.17 


1.00 








Q b h 2 


-0.67/0.00 


-0.37/0.01 


-0.05/0.01 


0.03/0.00 


0.01/0.00 


-0.87/0.02 


1.00/0.04 


1.00 






fi A 


-0.05/-0.02 


-0.06/-0.01 


-0.05/-0.01 


-0.01/0.00 


0.00/0.00 


-0.07/0.00 


0.05/0.01 


0.04/0.00 


1.00 




n s 


0.56/-0.01 


0.38/-0.04 


0.14/-0.03 


0.01/-0.01 


0.01/0.00 


0.85/-0.05 


-0.86/-0.20 


-0.83/0.03 


-0.13/0.00 


1.00 


A 


-/0.00 


-/0.00 


-/0.00 


-/0.00 


-/0.00 


-/0.00 


-/0.00 


-/0.00 


-/0.00 


-/0.00 1.00 



Note: We are not showing the foreground parameters. 



14 



15 - 
10 



5- 





<- -5 




-10 
-15- 
-20- 



16 18 20 22 24 

Z 

FIG. 7: Errors on the 7 parameter using SKAb(green) and 
SKAb+Planck (red, shorter error bars). In our analysis we 
assumed 7 was constant within bins of width Az — 2. The 
7 parameter is completely degenerate with xh and Sh, so 
in fact, what we can measure is "/xhSh (the values should 
then be multiplied by the same factor). Using Planck the 
degeneracy with the amplitude 8h is broken and we measure 
"fXH (smaller error bars). 
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FIG. 8: Errors on the j3 parameter using SKAb(green) and 
SKAb+Planck (red, shorter error bars). In our analysis we 
assumed (3 was constant within bins of width Az — 2. 



z = 7.5 and z = 9.5, with observations from 135 MHz to 
167 MHz over a frequency range of 32 MHz. 

In order to limit the number of parameters to be ex- 
tracted from the data, we parameterize the neutral frac- 
tion into two bins of width Az = 1.0. We consider this 
case as a direct comparison with previous estimates on 
the literature from Ref. (l3| . Note, however, that the as- 
sumed value for xh is possibly too high at the redshifts 
considered, if the optical depth to reionization, r is about 
0.1 as indicated by the WMAP3 analysis [13|. 

As we go to such low redshifts, the sky temperature 
decreases, so that experimental noise is lower (at z = 15 
T sky ~ 1500X while at z = 8 T sky ~ 340if). This helps 
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FIG. 9: Errors on the a parameter using SKAb(green) and 
SKAb+Planck (red, shorter error bars). In our analysis we 
assumed a was constant within bins of width Az — 2. 



first generation instruments such as LOFAR and MWA 
to possibly target these redshift ranges for observations 
and to make reasonable estimates on the ionization frac- 
tion, once cosmological parameters are fixed with Planck 
information; without some priors on the cosmological pa- 
rameters, neither LOFAR nor MWA will make mean- 
ingful measurements of the ionization fraction since it 
is degenerate with the normalization of the power spec- 
trum. With CMB information, all experiments, even the 
first generation ones, will establish the ionization fraction 
with a reasonable accuracy in several redshift bins. 

There are, however, certain complications and assump- 
tions. While the noise temperature is low, for a given 
interferometric array, however, low redshift observations 
are complicated by a corresponding increase in the beam 
size. This leads to a less number of modes at a given 
multipole to make anisotropy power spectrum measure- 
ment; this can be described as an increase in the cos- 
mic variance. Moreover, the boost provided by 7 when 
the signal is observed in absorption is no longer available 
since 7 = 1 with Tk >> T 7 and the amplitude could be 
smaller if xh is less than one, as the Universe is expected 
to begin, or even end the reionization process within this 
redshift range. Due to these reasons, the improvement in 
the parameter constraints at low redshifts with a given 
21-cm array is not as large as one might expect at first. 

In this limited scenario, in addition to constraints on 
the astrophysical parameters, the combination of 21- 
cm interferometric data with Planck does lead to cer- 
tain improvements on the cosmological parameters; ex- 
periments such as SKA and MWA5000 can improve 
cosmological parameters above the precision of Planck 
alone. These parameters include n s , tilt of the primor- 
dial power spectrum, that is improved by ~ 10% with 
MWA5000 and 40% with SKA. To the extent we can 
compare with Refs. [T2T [l3l|. we find generally consistent 
results, though we do expect differences as the two de- 
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FIG. 10: Marginalized elliptical error regions for pairs of model parameters from table ["VTl involving yiXH(zi) and either Qi,h 2 
or Qmh 2 . The contours correspond to a 68% likelihood on the joint parameters, with others marginalized over, using SKAb 
(dotted) and SKAb+Planck (solid). 



scriptions are different in the treatment of foregrounds. 
While we have utilized an explicit model for the fore- 
grounds, with their parameter uncertainties marginalized 
over when considering astrophysical and cosmological pa- 
rameter measurements, Ref. [12j ] assumed a foreground 
powerspectrum of the shape k~ 2 , with a free normaliza- 
tion, while Ref. [l3| removed foreground power at each 
three-dimensional Fourier mode under the assumption 
that the power in foregrounds is significantly less than 
the power in the 21-cm fluctuation and by correcting 
for the instrumental sensitivity. It is unclear which of 
these methods, including the one we use here where we 
marginalize over parameters related to the foreground 
model following Ref. [lj| , is the appropriate approach to 
reduce foreground contamination when analyzing 21-cm 
data. In any case, given the large amplitudes expected for 
the foregrounds, it is clear that any strong assumptions 
related to the foreground removal can have an impact on 



forecasting the ability to measure parameters with low 
frequency radio interferometers. 



C. When Tk >> T 1 and fluctuations in the 
ionization fraction are included. 

At the low redshifts considered in the previous analysis, 
reionization should already be well under way. Therefore, 
we now consider a fiducial model for the same frequency 
range as above, when more than 50% of the Universe 
is already ionized and the fluctuations in the ionization 
fraction are large. Table I Villi show the parameter con- 
strains in this case. Note that the size assumed for the 
bubbles is quite large as we followed the model in Ref. [|[ • 
This means that most of the signal from the perturba- 
tions in the ionization fraction will be smoothed out for 
I > 550 at z ~ 8 and I > 4700 at z ~ 9. Therefore, 
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TABLE VII: Forecasted l-cr uncertainties when Ts » T-, and b XH ~ (no significant fluctuations in the neutral fraction). 





XH 1 


XH2 


Q, m h 2 


n b h 2 




71 s 


Sh x 10 s 


Values 


0.95 


1.0 


0.127 


0.0223 


0.76 


0.951 


6.229 


SKAb 


0.30 


0.32 


0.030 


0.0122 


0.001 


0.008 




SKA 


0.51 


0.54 


0.051 


0.0208 


0.003 


0.014 




MWA5000 


0.62 


0.65 


0.066 


0.0255 


0.014 


0.027 




LOFAR 


30.1 


31.0 


3.10 


1.12 


0.35 


0.825 




MWA 


10.8 


11.3 


1.38 


0.48 


0.62 


0.919 




Planck 






0.0023 


0.00017 


0.011 


0.0047 


0.03 


SKAb + Planck 


0.012 


0.013 


0.0017 


0.00017 0.0013 


0.0025 


0.03 


SKA + Planck 


0.019 


0.020 


0.0017 


0.00017 


0.0031 


0.0030 


0.03 


MWA5000 + Planck 


0.044 


0.047 


0.0017 


0.00017 0.0086 


0.0039 


0.03 


LOFAR + Planck 


0.063 


0.087 


0.0023 


0.00017 


0.011 


0.0047 


0.03 


MWA + Planck 


0.082 


0.160 


0.0023 


0.00017 


0.011 


0.0047 


0.03 



Note: we used a frequency interval between z = 7.5 and z = 9.5 or 135 MHz to 167 MHz (we allow for the full cosmological evolution 
within this interval and also considered the expected foregrounds at these frequencies with their uncertainties marginalized over in the 
Fisher matrix formalism), xjj\ is the value of the neutral fraction between 7.5 < x < 8.5, while Xjj2 corresponds to the interval 
8.5 < z < 9.5, with both taken to be a constant across these redshift bins. We assumed a total of 2000 hours of observation on two places 

in the sky and a resolution of 0.1 MHz. 
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FIG. 11: The 21-cm power spectrum at z=8 (blue solid line) 
and z=9 (red dashed line) when b XH ~ (values from table 

rvTB . 



although the large bias factor (b XH ) should boost the sig- 
nal, this effect will be erased on small scales. As shown 
in Figure [P2l the angular power spectrum has character- 
istics features that depend on the angular scale at which 
density fluctuations begin to dominate the fluctuations 
in the ionization fraction. 

The overall results on forecasted parameter uncertain- 
ties are generally worse than the previous case as can be 
seen from the Table IVllTl This reduction in the precision 
of parameter estimates is due to the increase in the num- 
ber of parameters in this model (four more than the case 
where fluctuations in the neutral fraction are ignored) 
and a decrease of the signal power on small scales. Even 



with an increase in parameters, 21-cm interferometers 
when combined with Planck can measure the ionization 
fraction. The precision is better with experiments such 
as MWA5000 and SKA as they not only measure xh and 
b XH , but also parameters such as R XH in several redshift 
bins. Since the perturbations in the ionization fraction 
dominate the signal on scales larger than the bubble size, 
it may be easier for experiments to measure the product 
%h x b XH as a single parameter. This is specially true 
for first-generation instruments such as MWA and LO- 
FAR since these experiments will have problems measur- 
ing the signal from density fluctuations on scales smaller 
than the bubble size due to limitations on the beam size 
at the angular scales of interest. For the single parameter 
x H xb XH , we find errors of 1.54 (1.53) for LOFAR (MWA) 
at 7.5 < z < 8.5 and 0.20 (0.26) at 8.5 < z < 9.5. Since 
LOFAR captures some information from small scales, it 
is more capable of establishing the neutral fraction xh 
independent of the bias factor. 

Although we are considering the regime when Tg >> 
T 7 it might be difficult to verify a priori if we are actually 
in the regime where one should only consider fluctuations 
in the gas density or the neutral fraction. Therefore, 
when doing the data analysis, we should allow for some 
freedom in the other astrophysical parameters, Rl v , ol, 
(3 and 7 (as defined in eq. |2"5"|) . Note however that 7 will 
be completely degenerate with xh while the ae~ k Rl v/ 2 
term is quite small at these redshifts (in fact it is zero in 
our model), so that the only other parameter we might 
be able to constrain is (3. If f3 can be shown to be unity 
from the data alone, then this will serve as an indication 
that Ts >> T 7 ; the same can be tested with a direct 
measurement of the mean 21-cm brightness temperature 
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TABLE VIII: Forecasted 1-cr uncertainties when Ts >> T 7 and b XH is large 





XjJl 


XH2 


X H 1 




R XH1 (Mpc) R x 


„ n (Mpc) 


n m h 2 


n b h 2 




77 s 


5 H x 10 5 


Values 


0.2 


0.4 


-14.0 


-5.7 


50 


6 


0.127 


0.0223 


0.76 


0.951 


6.229 


SKAb 


0.04 


0.08 


0.42 


0.04 


3.9 


0.4 


0.020 


0.007 


0.0025 


0.018 




SKA 


0.11 


0.23 


0.58 


0.11 


11.5 


1.3 


0.058 


0.022 


0.0048 


0.040 




MWA5000 


0.19 


0.40 


1.07 


0.65 


29.3 


3.5 


0.145 


0.047 


0.017 


0.174 




LOFAR 


8.2 


16.7 


35.2 


9.0 


936 


111 


4.5 


1.70 


0.30 


3.01 




MWA 


4.1 


8.6 


36.8 


28.0 


889. 


110 


4.4 


1.23 


0.70 


7.39 




Planck 














0.0023 


0.00017 


0.011 


0.0047 


0.03 


SKAb + Planck 


0.004 


0.009 


0.37 


0.04 


0.53 


0.04 


0.0019 


0.00017 


0.002 


0.0041 


0.03 


SKA + Planck 


0.006 


0.015 


0.50 


0.08 


0.71 


0.05 


0.0021 


0.00017 


0.004 


0.0045 


0.03 


MWA5000 + Planck 


0.011 


0.044 


0.71 


0.36 


1.12 


0.11 


0.0022 


0.00017 


0.009 


0.0046 


0.03 


LOFAR + Planck 


0.12 


0.32 


30.1 


3.7 


44.0 


1.31 


0.0023 


0.00017 


0.011 


0.0047 


0.03 


MWA + Planck 


0.32 


1.17 


22.4 


13.2 


23.3 


3.1 


0.0023 


0.00017 


0.011 


0.0047 


0.03 



Note: We considered two separate bins in redshift with parameters with index 1 corresponding to the redshift interval 7.5 < z < 8.5, 
while parameters with index 2 corresponding to 8.5 < z < 9.5. See note in Table IVlTl for specifications used in these calculations. 
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FIG. 12: The 21-cm power spectrum at z=8 (blue solid line) 
and z=9 (red dashed line) when fluctuations in the ionization 
fraction are included (values from table IV1TT)) . 



before (with 50 Mpc and 6 Mpc) for the same redshift 
intervals (following 16]). Note that both LOFAR and 
MWA can achieve much better constraints in xh and 
b XH , independent of each other, when compared to the 
case where bubble fluctuations are larger. While this 
may look surprising, it is not unexpected. With a small 
bubble size, the damping scale moves to a small angular 
scale or a large multipole. This leads to extra power at 
multipoles between a 10 3 and 10 4 ; This difference in the 
amplitude of C; curves between Figs. [TT] and [13] lead to 
additional information on the parameters with the errors 
improved by the overall increase in the signal-to-noise 
ratio. When combined with Planck, experiments such 
as SKA and MWA5000 make significantly better mea- 
surements of all parameters relative to the case where 
bubble sizes are larger, ft could be that these interfer- 
ometers will pin down parameters related to reionization 
at a certain epoch when bubble sizes have not grown sig- 
nificantly large. 



since when T s > T 7 , the 21-cm signal will be in the emis- 
sion; unfortunately, it is unlikely that the mean temper- 
ature will be established with anisotropy measurements 
alone to verify if Ts » T 7 due to X-ray heating during 
reionization. 

The possibility to test whether Ts >> T 7 from 21-cm 
data is shown in tableflXl While with LOFAR and MWA, 
combined with Planck, one cannot establish whether /?$ 
is one, 21-cm experiments with sensitivity at MWA5000 
and better will allow measurements of /3j parameters at 
the 20% level or better; with SKA, this is at the 5% 
level. Measuring /3$ and showing that it is unity would 
be a direct test of heating in the IGM beyond CMB, such 
as due to a background of X-ray photons. 

Finally, we consider in table IXl the situation when the 
bubble size is much smaller than the suggested values 



VI. SUMMARY 

Here, we have studied the prospects for extracting cos- 
mological and astrophysical parameters from the low ra- 
dio frequency 21-cm background due to the spin-flip tran- 
sition of neutral Hydrogen during and prior to the reion- 
ization of the Universe. We made use of the angular 
power spectrum of 21-cm anisotropies, which exists due 
to inhomogeneities in the neutral Hydrogen density field, 
the gas temperature field, the intensity of the Lyman-a 
radiation from first luminous sources that emit UV pho- 
tons, and the gas velocity. Instead of the usual simplified 
case where fluctuations in the 21-cm brightness temper- 
ature are considered to be only due to the gas density 
field, when the spin temperature of neutral Hydrogen is 
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TABLE IX: Forecasted 1-a uncertainties when Ts >> T 7 and b Xjl is large 







XH2 


X H 1 


Ux H 2 


R XH1 (Mpc) R x 


„„ (Mpc) 


Pi 


Pi 


n m h 2 


Q b h 2 


n A 


n 3 


5 H x 10 5 


Values 


0.2 


0.4 


-14.0 


-5.7 


50 


6 


1.0 


1.0 


0.76 


0.951 


6.229 






Planck 


















0.0023 


0.00017 


0.011 


0.0047 


0.03 


SKAb + Planck 


0.006 


0.026 


0.46 


0.26 


0.53 


0.04 


0.02 


0.05 


0.0019 


0.00017 


0.0045 


0.0041 


0.03 


SKA + Planck 


0.012 


0.051 


0.79 


0.52 


0.72 


0.05 


0.05 


0.11 


0.0021 


0.00017 


0.0082 


0.0045 


0.03 


MWA5000 + Planck 


0.051 


0.071 


3.43 


0.80 


1.13 


0.11 


0.29 


0.20 


0.0022 


0.00017 


0.0105 


0.0046 


0.03 


LOFAR + Planck 


1.1 


0.80 


83 


10 


44 


1.35 


6.5 


2.4 


0.0023 


0.00017 


0.011 


0.0047 


0.03 


MWA + Planck 


2.1 


1.37 


149 


21 


24 


5.65 


13.1 


8.6 


0.0023 


0.00017 


0.011 


0.0047 


0.03 



Same as table IVTlTl but also considering as an extra parameter to be extracted from the data (with the fiducial value of 1, by definition, 

since T s > T 7 ). 



TABLE X: Forecasted 1-a uncertainties when Ts » T 7 , b XH is large and R XH is small 







XH2 


bx H1 bx H2 


R XH1 (Mpc) R x 


„ 2 (Mpc) 




n b h 2 


H A 


n s 


S H x 10 5 


Values 


0.2 


0.4 


-14.0 -5.7 


1.0 


0.5 


0.127 


0.0223 


0.76 


0.951 


6.229 


SKAb 


0.01 


0.03 


0.012 0.006 


0.005 


0.014 


0.004 


0.0022 


0.0006 


0.003 




SKA 


0.03 


0.07 


0.018 0.011 


0.012 


0.038 


0.011 


0.0056 


0.0012 


0.005 




MWA5000 


0.06 


0.12 


0.038 0.027 


0.059 


0.238 


0.027 


0.0111 


0.0027 


0.010 




LOFAR 


1.20 


2.35 


0.70 0.48 


1.14 


5.18 


0.48 


0.15 


0.048 


0.17 




MWA 


0.66 


1.32 


1.61 1.30 


3.38 


12.60 


0.44 


0.15 


0.074 


0.36 




Planck 












0.0023 


0.00017 


0.011 


0.0047 


0.03 


SKAb + Planck 


0.002 


0.004 


0.012 0.006 


0.005 


0.014 


0.0015 


0.00017 


0.0006 


0.0020 


0.03 


SKA + Planck 


0.003 


0.006 


0.017 0.012 


0.011 


0.039 


0.0016 


0.00017 


0.0012 


0.0028 


0.03 


MWA5000 + Planck 


0.004 


0.008 


0.027 0.026 


0.052 


0.236 


0.0015 


0.00017 


0.0025 


0.0036 


0.03 


LOFAR + Planck 


0.011 


0.044 


0.30 0.47 


0.91 


5.1 


0.0023 


0.00017 


0.0106 


0.0047 


0.03 


MWA + Planck 


0.015 


0.110 


0.69 1.30 


2.03 


12.3 


0.0023 


0.00017 


0.0108 


0.0047 


0.03 



Same as table I Villi but assuming bubble size is smaller so that the fluctuations in the ionization fraction are basically dominating the 

signal on all scales. 



significantly larger than that of the CMB, and routinely- 
considered in the literature [U, [l6[ , we considered a gen- 
eral case where fluctuations are induced by a variety of 
sources during the era of reionization. The model in- 
cludes sources of 21-cm brightness fluctuations during an 
era when Lyman-a photons are present and due to vary- 
ing levels of neutral-fraction fluctuations. 

Using a Fisher analysis for a variety of upcoming and 
planned low-frequency 21-cm interferometers, we fore- 
casted uncertainties in parameters that describe both the 
underlying mass power spectrum and the global cosmol- 
ogy, as well as a set of simplified astrophysical parameters 
that connect fluctuations in the dark matter to those that 
govern 21-cm fluctuations. In addition to detector noise, 
we also marginalized over a model for the foregrounds at 
low radio frequencies; the limiting factor for 21-cm ob- 
servations are generally the residual noise level from fore- 
grounds. Using our general description for 21-cm bright- 
ness temperature anisotropics, we find large degeneracies 
between the cosmological parameters and the astrophys- 
ical parameters; such degeneracies are present to a less 



extent when strong assumptions are made with respect 
to the 21-cm spin temperature relative to the CMB tem- 
perature such as the limiting scenario that Ts >> T~ f \ 
but this decrease is primarily due to a reduction in the 
number of parameters to be extracted from the data and 
these parameters are all astrophysical quantities that are 
strongly degenerate with certain cosmological parame- 
ters. 

We have shown how the parameter degeneracies are 
broken when 21-cm measurements are combined with in- 
formation from the CMB, such as using anisotropy and 
polarization data expected from Planck. For upcoming 
low frequency radio interferometers, the overall improve- 
ment on the cosmological parameter estimates when com- 
bined with Planck is not significant. Interferometers such 
as LOFAR and MWA, however, will be able to mea- 
sure astrophysical parameters such as the neutral frac- 
tion at the tens of percent level. Experiments such as 
MWA5000 and SKA will improve cosmological parame- 
ter constraints beyond that of the Planck alone, while at 
the same time improving measurements of astrophysical 
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10 1 10 2 10 3 10 4 



quantities during the reionization era. These interferom- 
eters will also establish the reionization history, the onset 
of gas heating (by showing $ = 1 at some epoch) and 
also establish, for example, if Tk < T-y. 



FIG. 13: The 21-cm power spectrum (solid lines) and corre- 
sponding errors for LOFAR (dotted lines) and MWA (dashed 
lines) when fluctuations in the ionization fraction are included 
and R XfI is small (values from table [X). Thin red lines cor- 
respond to 2 ~ 9 while thick blue lines to z ~ 8. Note that 
the errors shown here take fully into account the number of 
modes available at each I. Note also that this is just for one 
frequency bin (we have a total of 320 bins that can be used 
to improve the measurements for the two redshift bins con- 
sidered) . 
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